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1. 


Introduction 


Satellite  communication  links  at  UHJF  can  be  subject  to  the  effects  of 
ionospheric  scintillations.  These  scintillations  cause  both  enhancements  and 
fading  about  the  median  level  as  the  radio  signal  transmits  the  disturbed 
ionospheric  region.  When  scintillations  occur  which  exceed  the  fade  margin, 
performance  of  the  communications  link  will  be  degraded.  Of  major  Importance 
is  the  estimation  of  the  occurrences  of  these  scintillations  that  result  in 
degradation  of  the  communications  link.  One  approach  to  this  estimation 
problem  consists  of  determining  the  probability  distribution  of  scintillations 
and  then  using  the  properties  and  parameters  of  this  distribution  to  obtain 
the  desired  estimates  related  to  fading.  The  Nakagami-m  distribution 
(Nakagami,  1960)  has  been  shown  to  be  a  useful  distribution  for  describing 
the  effects  of  scintillations  (Whitney,  Aarons,  Allen  and  Seeman,  1972). 

This  paper  discusses  the  statistical  procedures  that  are  applicable  in 
attempting  to  judge  the  goodness-of-f it ,  or  appropriateness  of  the  Nakagami-m 
distribution  to  a  data  set.  That  is,  it  discusses  procedures  for  determining 
if  a  given  data  set  can  be  considered  a  sample  of  data  generated  from  a 
Nakagami-m  distribution. 

While  the  underlying  problem  which  this  paper  addresses  did  arise 
from  an  investigation  of  scintillation  data,  the  statistical  techniques  are 
not  specific  to  this  problem.  The  techniques  are  applicable  to  any  data  set 
arising  as  either  independent  observations  or  as  a  stationary  time  series 
which  might  be  from  a  Nakagami-m  distribution. 


\  iy . .. 
j  p is ‘ • 
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2. 


Mathematical  Properties  of  the  Nakagami-m  Distribution 

The  probability  density  for  the  Nakagarai-m  distribution  is  normally 
given  as  an  amplitude  or  power  probability  density  function 


m 


m 


f  (s)  =  _ 

T  (m)_m 


m-1  ,  ms  . 

s  exp  (-  -g-  ) 


(2.1) 


where 

s  =  signal  power  (watts), 
fl  =  average  power, 

1/2  <  m  £  co 

and 

T(m)  =  gamma  function  of  m. 

For  modelling  scintillation  data  ra=l  in  (2.1)  is  referred  to  as  Rayleigh 
fading.  In  such  a  situation  the  probability  density  is 

fS(s)  =  §  exp  (2.2) 

This  distribution  (2.2)  is  usually  called  the  exponential  or  negative 
exponential  distribution.  (2.2)  is  related  to  the  classical  Rayleigh  distri¬ 
bution  when  we  consider  the  transformation  of  variables 

R2  =  S  (2.3) 


where  R  is  intensity.  The  probability  density  of  R  is  then 


fR(r)  -  -¥  exp  (-  -f- } 


(2.4) 


which  Is  the  classical  Rayleigh  distribution.  Also  for  scintillation  data 
m  1  refers  to  fading  more  severe  then  Rayleigh  fading. 
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The  nth  moment  of  the  random  variable  S  of  (2.1)  about  the  origin  is 


ES 


n 


fs(s)  ds  =  1  (-rcm)-^  > 


- ]  (2.5) 

(^)  (n+m) 


Fcr  scintillation  data  an  important  parameter  is  the  coefficient  of  variation 
or,  in  scintillation  jargon,  the  index.  Here 

1/2 

O  1 


=  (ES2  -  (ES)2) 

4  ES  U  /  m 


(2.6) 


In  (2.6)  we  use  u  and  O  to  represent,  respectively,  the  mean  and  standard 
deviation,  m  is  the  parameter  of  the  Nakagami-m  distribution  given  in  (2.1). 

3.  Relation  to  the  Gamma  Distribution 


The  Nakagami-m  distribution  as  given  in  (2.1)  is  related  to  the  more 
standard  gamma  distribution  whose  density  is  given  by 


fs(s) 


a- 1 

s 


r(  oox3 


exp 


(3.1) 


Notice  if  we  set  in  (3.1) 


a  =  m  and  X  =  Q/m 


(3.2) 


the  gamma  density  of  (3.1)  is  equal  to  the  Nakagami-m  density  of  (2.1). 

This  relationship  is  important  for  it  allows  us  to  use  the  extensive  theory 
developed  for  the  gamma  distribution  to  solve  problems  dealing  with  the 
Nakagami-m  distribution. 


6 


4.1  Ecdf  and  Kolmogorov-Smirnov  Test 


Say  X^,...,Xn  represent  a  sample  of  size  n.  Further,  say  we  wish  to  evaluate 
whether  this  sample  came  from  a  Nakagami-m  distribution.  To  begin  the  analy¬ 
sis  one  should  first  compute  the  empirical  cumulative  distribution  function 
(ecdf)  defined  for  arbitrary  x  as 


Fn  (x) 


#(X±  £  x) 
n 


(4.1) 


Figure  la  contain  an  ecdf  plot  of  10  random  observations  drawn  from  a  nega¬ 
tive  exponential  distribution  with  mean  5  (i.e.,  Rayleigh  fading,  m=l  and 
fl=  5  in  (2.1)).  The  10  observations  are:  8.15,  4.69,  2.17,  0.37,  16.69, 
0.06,  6.48,  2.63,  0.44,  0.89.  The  sample  arranged  in  order  of  magnitude  and 
with  the  corresponding  ecdf  values  are: 


Ordered  Observation 
Number  (i) 

Ordered 

Observation 

Ecdf 

Fn  (X) 

1 

0.06 

.10 

2 

0.37 

.20 

3 

0.44 

.30 

4 

0.89 

.40 

5 

2.17 

.50 

6 

2.63 

.60 

7 

4.69 

.70 

8 

6.48 

.80 

9 

8.15 

.90 

10 

16.69 

1.00 

The  next  step  is  to  plot  the  cumulative  distribution  function  (cdf)  for  the 
hypothesized  distribution  on  the  same  graph  with  the  sample  ecdf  and  then 
judge  if  the  ecdf  differs  significantly  from  the  hypothesized  cdf.  For  a 
continuous  random  variable  X  with  probability  density  F(x)  the  cdf  is  defined 


(A. 2) 


x 

F(x)  *=  J  f (y)  dy 

-  tO 

Say  for  the  present  example  the  hypothesized  distribution  is  the  exponential 
given  in  (2.2).  The  edf  for  this  distribution  is 


F(x)  -  1  -  e  x/Q 

and 

F(x)  =  0 


for  x  >  0 


for  X  <  0 


(4.3) 


Recall  this  distribution  represents  Rayleigh  fading  (i.e,  m=l  in  (2.1)). 

Two  situations  present  themselves  here.  First,  the  values  of  all  the  parame¬ 
ters  of  the  hypothesized  distribution  are  known.  For  our  example  the  only 
parameter  is  ft.  Figure  lb  contains,  in  addition  to  the  ecdf  of  the  observa¬ 
tions,  the  cdf  of  (4.3)  with  ft=5.  The  second  situation  is  when  the  values  of 
some  of  the  parameters  are  not  known.  Figure  1c  contains,  in  addition  to  the 

ecdf  of  the  ten  observations,  the  cdf  of  (4.3)  where  ft  is  replaced  by  an 

A 

estimate,  ft,  of  it  which  is  the  sample  mean  X  =  4.26.  In  general  the  unknown 
parameters  should  be  replaced  with  efficient  estimates  -  e.g.,  minimum  vari¬ 
ance  estimators  or  maximum  likelihood  estimators.  However,  estimators  ob¬ 
tained  by  the  method  of  moments  are  also  often  used  for  "quick  computations". 
For  our  example  the  moment  estimator  is  also  the  minimum  variance  and  maximum 
likelihood  estimator. 

To  Judge  the  significance  of  the  difference  between  the  cdf  and  the  ecdf  the 
investigator  has  two  possibilities.  The  first  is  simply  to  judge  informally 
if  the  difference  is  too  large.  For  example,  the  investigator  can  compute 

F  (x)  -  F  (x)  (4.4) 

n 

for  a  variety  of  x's  and  make  a  judgement  concerning  their  magnitudes.  In 
this  situation  the  investigator  is  usually  asking  the  question  "Are  the 
differences  in  (4.4)  of  any  practical  significance?  The  second  procedure 
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consists  of  using  a  formal  statistics  test  of  significance  -  viz.,  the 
Kolmogorov-Smirnov  test  (see  Dixon  and  Massey,  1969,  p.  345).  This  test  con¬ 
sists  of  computing 

D  =  sup  If  (x)  -  F  |  (4.5) 

X  x ' 

and  rejecting  the  hypothesized  distribution  as  the  true  distribution  if  D  of 
(4.5)  exceeds  a  critical  value,  say  d  The  value  d^  is  selected  to  produce 
a  test  of  level  of  significance  equal  to  a-  i.e.,  it  is  selected  so  that 
there  is  an  a  chance  of  D£dQ  if  the  hypoth4sized  distribution  is  the  true 
distribution.  Alternatively  this  test  consists  of  adding  d  to  all  values  of 
F^(x).  The  results  of  such  a  computation  are  shown  in  Figure  Id.  (Note  in 
Figure  Id  F^  (x)  —  d  is  forced  to  lie  between  0  and  1.  Also  for  this  figure 
dQ  =  .41  for  n  =  10  and  a=  .05).  If  any  of  the  cdf  is  outside  the  band,  the 
hypothesized  distribution  is  rejected  as  the  true  distribution  at  the  o  level 
of  significance. 

This  last  version  of  the  Kolmogorov-Smirnov  test  can  also  be  used  to 
produce  a  confidence  interval  or  region  for  the  underlying  distributions  cdf. 
Any  cdf  lying  completely  in  the  band  F^  (x)  —  dQ  is  an  acceptable  cdf  at  the 
100  (1-c)  percent  level  of  confidence.  For  example,  the  shaded  area  in  figure 
Id  consists  of  the  95  percent  confidence  region  for  the  10  random  observations 
given  above. 

The  values  dQ  of  the  Kolmogorov-Smirnov  test  depend  upon  the  de¬ 
sired  level  of  significance  (or  desired  confidence  level)  and  the  sample  size. 
One  table  of  d  values  is  given  in  Dixon  and  Massey  (1969).  When  the  sample 
size  of  n  independent  observations  exceeds  30  the  following  values  of  dQmay 
be  used: 

Significance  Level  Confidence  Level  da 


.10 

.90 

1.22/lS  n 

.05 

.95 

1.36/Vn 

.01 

.99 

1.63/i/  n 
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4 • 2  Special  Considerations  for  Application  nj  Kolmogorov-Smirnov  Test 


The  Kolmogorov-Smirnov  test  as  described  above  is  applicable  to  situations 
where  we  have  samples  consisting  of  independent  observations  and  all  the 
parameters  of  the  hypothesized  distribution  are  given  explicitly.  If  the 
observations  are  not  independent,  as  will  happen  when  we  have  a  time  series, 
then  we  have  two  possible  procedures.  First,  the  test  can  be  applied  to  only 
a  subset  of  the  observations.  One  way  of  obtaining  this  subset  is  to  compute 
the  autocorrelation  function,  find  the  period  or  lag  that  corresponds  to  a 
zero  autocorrelation  (say  it  is  period  k)  and  use  every  kt’u  observation  in  the 
Kiomogorov-Smimov  test.  Alternatively,  the  ecdf  of  (4.1)  can  be  computed 
using  all  the  observations,  but  the  values  of  d  should  be  multiplied  by  A- 
This  will  effectively  reduce  the  sample  to  n/k  independent  or  uncorrelated 
observations  without  the  loss  of  any  information  obtainable  from  the  full  set 
of  n  available  observations. 

If  any  of  the  parameters  are  not  known  then  they  must  first  be  estimated  before 
the  Kolmogorov-Smirnov  test  can  be  applied.  For  the  Nakagami-m  distribution 
of  (2.1)  the  two  parameters  that  need  to  be  estimated  are 

ft  and  m 


These  can  be  estimated  by  the  method  of  maximum  likelihood  or,  if  the  sample 
is  large,  there  should  be  little  loss  in  efficiency  if  the  parameters  are 
estimated  by  the  method  of  moments.  The  moment  estimates  are: 


and 


a  _ 

Q  =  x 


-  2 

A  X 

m  =  — 2 
s 


(4.6) 


(4.7) 


In  (4.5)  and  (4.6)  x  and  s  are  the  sample  mean  and  variances,  respectively, 
where 


-  Ex 
x  =  — 
n 


2  =  E (x-x) ' 
n-1 
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When  the  parameters  are  estimated  and  then  the  Kolmogorov-Smirnov  test  is 
applied,  the  resulting  test  is  conservative.  That  is,  the  true  level  of 
significance  is  smaller  than  the  nominal  or  stated  level. 


4 . 3  Probability  Plotting 


If  the  above  procedure  leads  to  rejection  of  the  hypothesis  that  the 
Nakagami-m  distribution  "fits"  the  data,  the  next  item  in  the  analysis  is  to 
determine  where  and  how  the  model  deviates  from  the  data.  Certain  deviations 
may  not  be  considered  to  be  of  practical  significance.  For  example,  it  may 
not  be  a  serious  lack  of  fit  if  the  data  deviates  from  the  model  only  for  the 
tail  observation  (say,  less  than  2nd  percentile  or  greater  than  98  percentile). 
However,  other  deviations  may  be  considered  very  serious  and  would  render  the 
Nakagami-m  model  useless.  it  is  important  that  the  investigator  knows  the 
statistically  significant  deviations  and  knows  if  they  are  of  practical  im¬ 
portance  . 


One  useful  way  of  determining  where  the  model  deviates  from  the  data  is  to  em¬ 
ploy  probability  plotting.  Probability  plotting  for  the  present  problem  is 
the  plotting  of  the  ordered  observations  from  a  sample  versus  the  inverse  of 
the  cdf  of  the  Nakagami-m  distribution  of  (2.1).  Usually  fi  is  set  equal  to 
unity  in  the  plotting.  Specifically,  say  x^  ^  .  .  .  Cx(n)  represents  the 
sample  ordered  from  the  smallest  to  largest.  Next,  say  the  Nakagami-m 
distribution  with  Q=1  has  cdf  (G(z)  equal  to 


z  m 

G  (z)  =  f  y11  exp  (-my)  dy 

Jo  7 


Further  say  we  redefine  the  ecdf  to  be 

F(x  )  - 
U(k)'  n 


(4.8) 


(4.9) 


The  probability  plotting  is  a  plot  on  1 inear-by-1 inear  paper  of 


X(i)  on  G  1(F  (x(i)^ 


(4.10) 
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If  the  Nakagami-m  distribution  is  the  "true"  distribution  the  plot  given  by 
(4.10)  is,  within  sampling  fluctuations,  a  straight  line  through  the  origin 
with  slope  equal  to  Q.  Deviations  from  a  straight  line  indicate  where  the 
deviation  from  the  model  exists. 

Notice  in  order  to  implement  the  above  m  must  be  known. 

One  possibility  for  determining  m  if  it  is  not  known  a  priori  is  to  use  the 
estimator  m  given  by  (4.7).  Also  note  that  the  inverse  of  the  cdf  defined  in 
(4.8)  cannot  be  given  in  closed  form.  However,  Wilk,  Gnanadesikan  and  Huyett 
(1962)  supply  tables  and  an  outline  for  a  computer  program  which  can  be  used 
to  obtain  these  inverses. 
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5.  Chi-Square  Goodness-of-Fit  Test 

We  recommended  the  Graphical  Analysis  coupled  with  the  Kolmogorov  Smirnov  test 
described  above  as  the  preferred  technique  for  determining  the  appropriate¬ 
ness  of  the  Nakagami-m  distribution  to  a  set  of  data.  However,  here  may  be 
situations  where  the  Kolmogorov-Smirnov  test  may  not  be  applicable.  For  ex¬ 
ample,  the  researcher  may  know  a  priori  that  the  extreme  tails  of  the  data 
(below  2nd  percentile  and  above  98th  percentile)  will  not  be  well  approximated 
by  the  Nakagami-m  distribution.  This  could  be  due  to  accuracy  limitations  of 
the  measurement  instrument.  In  such  a  situation  the  investigator  may  want  to 
censor  (i.e.,  remove)  the  tails  of  the  data  and  not  enter  these  into  a  formal 
statistical  inference  test.  The  ecdf  defined  in  (4.1),  the  plot  of  the  ecdf 
(such  as  in  figure  la) ,  and  the  probability  plotting  as  described  in 
section  4.3  are  still  valid  and  useful.  However,  the  Kolmogorov-Smirnov  test 
is  not  valid  on  censored  data.  The  chi-square  goodness-of-f it  test  is  appro¬ 
priate  in  this  situation  as  a  statistical  inference  test.  We  suggest  the  test 
should  be  performed  as  follows: 

(1)  Decide  upon  appropriate  values  of  m  and  £3.  These  may 
be  known  as  apriori  or  estimated  from  the  data.  The  appro¬ 
priate  method  for  obtaining  these  is  by  use  of  the  method 

of  maximum  likelihood  on  the  censored  data.  However,  if  the 
sample  is  large  the  method  of  moment  estimates  (see  (4.6)) 
and  (4.7))  should  be  sufficient. 

(2)  Using  the  values  of  m  and  (1  obtained  from  (1)  find  the 
values  Sg, . . .S^ , . . .S^q  which  divide  the  distribution  (2.1) 
into  20  equal  probability  sections.  Note  Sq=0,  S1  is  de¬ 
termined  such  that 

.05  =  fl  fs  (s)  ds 
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S0  is  determined  such  that 
S2 

.10  “  fg  (s)  ds^  . . . 

0 

and  - 'J* .  This  step  produces  20  intervals  or  categories 
each  with  probability  .05. 

(3)  Compute  the  expected  values  for  each  of  the  20  categories. 
These  will  all  equal  n(.05). 

(4)  Classify  each  observation  into  one  of  the  twenty  categories. 
The  frequencies  in  these  20  categories  can  be  represented  by 

f^,  ....  where  n  =  Ef^. 


(5)  Compute  the  chi  square  statistic 

X2  =  E(f  -  . 05n) 2 / ( . 05n)  (5.1) 

2 

(6)  Compare  the  X  value  of  (5.1)  with  the  appropriate  critical 
chi  square  value  obtained  from  the  chi  square  distribution  with 
19  degrees  of  freedom  if  m  and  £1  were  estimated,  or  obtained  from 
the  chi  square  distribution  with  17  degrees  of  freedom  if  both 

m  and  Q  were  estimated. 

If  the  data  contains  dependent  observations  (as  in  a  time 
2 

series)  then  the  X  value  of  (5.1)  should  be  divided  by  k,  i.e.. 


X  =  X  /  k 
new 


where  k  is  the  period  or  lag  corresponding  to  a  zero  correlation 

2 

in  the  oroginal  data  (see  section  4.2).  The  X  value  is  now 

new 

compared  to  the  critical  chi  square  values  obtained  from  the  chi 
square  tables. 
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6.  Test  for  Changing  m  Values 

One  problem  which  the  authors  have  had  to  address  when  attempting  to  evaluate 
the  appropriateness  of  the  Nakagami-m  to  time  series  is  the  problem  of  chang¬ 
ing  m  values.  That  is,  while  the  Nakagami-m  distribution  may  be  an  appropriate 
model  for  the  data,  the  actual  value  of  m  is  not  constant  over  the  entire  data 
set.  In  some  of  these  situations,  the  number  and  locations  of  the  segments 
that  have  different  m  values  may  be  known.  In  this  section,  we  present  a 
large  sample  test  which  can  be  used  to  test  the  equality  of  the  Nakagami-m 
values  for  t  segments  of  data. 

6. 1  Mathematical  Statement  of  the  Problem 

Say  we  have  t  sets  of  Independent  observations  each  from  a  Nakagami-m 
distribution.  The  m  values  for  the  segments  are 

^1  *  m2  *  *  *  *  * 

The  problem  is  to  test  for  the  equality  of  these  m  values.  That  is,  we  want 
to  test  the  hypothesis 

H:  nij  =  =  ...  *  m  ■  m.  (6.1) 

Equivalent  to  this  hypothesis  is  the  hypothesis 

»1;  SA1  =  S42  =  •”  =  S4t  =  S4  (6*2) 

where  S.,  for  i  =  1,  . . . ,  t  is  the  S.  value  defined  in  (2.6)  for  the  i^1 
si  4 

segment.  (Recall  S,  =  1/  / m.) .  The  test  we  present  in  the  following  tests 
directly  the  hypothesis  H  of  (6.2). 
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6.2  Notation  -  Statistical  Resulta 


Say  we  have  independent  observations  from  segment  i  for  1=1,  ....  t.  From 
each  segment  we  compute  the  sample  mean,  sample  standard  deviation,  and  sample 
estimate  of  S^.  These  are,  respectively, 

Xl  «  E  (Xlj)/n1  (6.3) 


n-1 


(6.4) 


and 


(6.5) 


Here  represents  the  jth  observation  for  the  1th  sample,  j-1,  ...,  n^, 
1*1,  ....  t.  For  large  samples 


(6.6) 


where  E  represents  the  expected  value  operator.  Further  for  large  samples 

a 

the  standard  error  of  for  i*l,  ...»  t  is 


'41 


1 

1 

U41  -“21 

+  U21 

-  u3i 

ni 

— 

mi 

4W21 

U2i  Ui_ 

1/2 


(6.7) 


where 


Vl=EXij 


(6.8) 


17 


and 


U^i  -  E(X  -  it±)  *  for  l  y-  0 


(6.9) 


The  sample  estimate  of  the  standard  errof  of  is 


4i 


Here 


1 

ni 

A  2 

(S4i> 

A  ^  A  A 

U4l  "  U2i  +  ^2i  -  U3i 

4u2i  2  ^2  A. 

_ 

U  .  = 

1 


E(xn  -  V 

"i 


for  i  =  1,  ....  t. 


1/2 


(6.10) 


Further  for  large  samples  the  S  are  approximately  normally  distributed. 
See  Rao  (1973,  Chapter  6)  for  proofs  of  the  above  assertions. 

If  the  m^  values  are  all  equal  (that  is,  hypothese  H  of  (6.1)  and  of 

(6.2)  are  correct,  then  an  estimate  of  the  common  value  m  is 

A  n  2 
E(n  S  ag  ) 

*  1  41  S4i 


E(n./o  *2  ) 

1  b4i 


(6.11) 


6.3  The  Test 


Given  that  H  :  S 


41 


.  *  S.  is  correct  then 
4 1 


E.- 

l 


A  A 

(S4U  -  V 

0*2 

b4i 


(6.12) 
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i 


m 


is  approximately  distributed  as  a  chi  square  variable  with  t-1  degrees  of 
freedom  for  large  samples.  Rejection  of  H1  at  the  a  level  of  significance 
follows  if  the  statistic  of  (6.12)  exceeds  the  upper  o  value  of  the  chi- 
square  distribution  with  t-1  degrees  of  freedom  (Rao,  1973,  p.  389). 

6 . 4  Further  Comments 

In  addition  to  the  m  values  varying  from  segment  to  segment  the  ft  value  of 
the  Nakagami-m  distribution  may  also  vary.  This  value  is  the  mean  of  the 
distribution  so  an  appropriate  test  to  test  the  hypothesis 


is  the  analysis  of  variance  test.  Because  the  mean  of  the  Nakagami-m 
distribution  is  proportional  to  its  standard  deviation  the  analysis  of 
variance  on  the  logs  of  the  data  may  be  a  more  appropriate  analysis  than  an 
analysis  of  variance  of  the  original  data  (Dixon  and  Massey,  1969,  Chapter  16). 

As  was  stated  a  number  of  times  above,  the  test  for  equality  of  the  m  values 

assumes  independent  observations.  If  we  are  dealing  with  a  time  series  then 

the  sample  sizes  should  be  reduced  or  other  adjustments  should  be  made  to 

t  h 

reflect  this  (see  section  4.2).  One  possibility  is  to  use  every  k  obser- 

th  * 

vation  in  the  i  segment  for  i-1,  ...,  t  where  k  Is  the  period  or  lag 

th  ^ 

corresponding  to  zero  autocorrelation  for  the  i  segment.  Then  only  n^/k^ 
observations  will  be  used  in  each  segment.  Alternatively,  all  n^  observations 
can  be  used  by  n.  in  formulas  (6.7)  and  (6.10)  should  be  replaced  by  n,/k.. 

1  K  1 
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